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Ras GTPases are lipid- anchored G proteins which play a fundamental role in cell signaling pro- 
cesses. Electron micrographs of immunogold-Iabeled Ras have shown that membrane-bound Ras 
molecules segregate into nanocluster domains. Several models have been developed in attempts 
to obtain quantitative descriptions of nanocluster formation, but all have relied on assumptions 
such as a constant, expression-level independent ratio of Ras in clusters to Ras monomers (clus- 
ter/monomer ratio). However, this assumption is inconsistent with the law of mass action. Here, 
we present a biophysical model of Ras clustering based on short-range attraction and long-range 
repulsion between Ras molecules in the membrane. To test this model, we performed Monte Carlo 
simulations and compared statistical clustering properties with experimental data. We find that we 
can recover the experimentally-observed clustering across a range of Ras expression levels, without 
assuming a constant cluster/monomer ratio or the existence of lipid rafts. In addition, our model 
makes predictions about the signaling properties of Ras nanoclusters in support of the idea that 
Ras nanoclusters act as an analog-digital-analog converter for high fidelity signaling. 



I. INTRODUCTION 

Plasma membrane heterogeneity is a key concept in 
molecular cell biology due to its role in protein sorting 
and specificity of signaling [J, [3, 0] ■ Although the diver- 
sity of the membrane's lipid components is partly respon- 
sible for this heterogeneity [1], the role played by mem- 
brane proteins is less well understood. Members of the 
Ras protein superfamily Q have been observed to form 
dynamic, non-overlapping domains called nanoclusters in 
the inner leaflet of the plasma membrane 0, H, H, fiot . 
While the lateral segregation of Ras may provide evi- 
dence towards the existence of small, dynamic rafts [ll[ , 
the definition and even existence of rafts remains dis- 
puted [TH ■ In addition to its connection to the lipid-raft 
concept, Ras has attracted immense interest due to its 
fundamental role in a multitude of cellular processes, in- 
cluding cell proliferation, survival, and motility. Most 
importantly, Ras genes are found to be mutated in 30% 
of human cancers |l3l . [l^ . [Tsj , making their products ex- 
tremely important therapeutic targets [lB| . While the in- 
tracellular biochemistry of Ras genes is well documented, 
the biophysical mechanism and role of Ras clustering in 
the plasma membrane remains little understood. 

Ras GTPases are small (21 kDa), lipid-anchored pe- 
ripheral membrane proteins involved in signal transduc- 
tion [31 . Three Ras isoforms H-Ras, K-Ras and N-Ras 
are expressed in all mammalian cells. These isoforms 



contain a conserved G-domain which binds guanine nu- 
cleotides [13. Ras effectively acts as a molecular switch 
for the signal, with "on" (GTP-bound) and "off" (GDP- 
bound) states, the former promoting an association with 
and activation of effector proteins. Although nearly iden- 
tical with respect to their catalytic and effector-binding 
properties, H-Ras, N-Ras and K-Ras have very different 
biological roles. This functional distinction is believed 
to result at least in part from the differential membrane 
compartmentalization of Ras isoforms [l^, [l^ . The dif- 
ferent distribution of Ras proteins in cellular membranes 
dictates unique spatio-temporal patterns of activation 
of effector pathways. A classical example of a path- 
way involving Ras is the Ras-Raf-MEK-ERK pathway, 
a mitogen-activated protein kinase (MAPK) cascade in- 
volved in cell proliferation, differentiation, and apoptosis. 
In this pathway, the epidermal growth factor receptor 
(EGFR), a receptor tyrosine kinase, is stimulated. This 
leads to recruitment and activation of guanine nucleotide 
exchange factors (GEFs) which, by interacting with the 
Ras G-domains, promote the exchange of GDP for GTP 
(irt and lead to Ras activation. Ras-GTP activates pro- 
tein kinase Raf and initiates the phosphorylation cas- 
cade, ultimately leading to double phosphorylated ERK 
(ERKpp) which then travels into the nucleus and phos- 
phorylates transcription factors (20| . Among other pur- 
poses, such cascades can lead to a massive amplification 
of the original signal [2^ . 
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FIG. 1: Experimental immunoEM data and statistical clus- 
tering analysis. (A) Electron micrograph of immunogold- 
labeled Ras domain (GFP-tH where tH is minimal plasma 
membrane targeting motifs of H-Ras) in an in vitro plasma 
membrane sheet. Scale bar is 100 nm. (B) Correspond- 
ing point-pattern analysis (red) and 99% confidence interval 
(black). ©Prior et al. (2003), originally published in The 
Journal of Cell Biology. doi:10.1083/jcb. 200209091 0. 



Experimental evidence for the formation of nanoclus- 
ters (termed clusters from now on) is provided by in vivo 
and in vitro experiments. Fluorescence resonance energy 
transfer (FRET) studies show that activation by EGF 
leads to significant decrease in Ras lateral diffusion, sug- 
gesting the existence of Ras-GTP clusters [2l[. A very 
similar result was obtained by single-molecule fluores- 
cence microscopy, where GTP-binding of Ras leads to 
slowly diffusing active Ras molecules [23 ■ Single particle 
tracking (SPT) studies of fluorescently labeled Ras have 
also demonstrated transient immobility of Ras (lasting 
less than Is) with high temporal resolution, interspersed 
with periods of free Brownian motion [23j . Furthermore, 
spatial statistics of fluorescently labeled Raf have shown 
that Ras and Raf cluster together [23|. It is therefore 
believed that active Ras forms signaling platforms which 
recruit and activate Raf. As signaling platforms are Ras- 
isoform specific, the signal diversity observed between H- 
Ras, K-Ras and N-Ras is in part the result of differential 
clustering properties in these isoforms 01. 

Direct evidence for protein clustering in a membrane 
can be obtained from high-resolution electron microscopy 
(EM). However, Ras is too small and not electron dense 
enough to be observed directly. To circumvent this prob- 
lem, Prior et al. used GFP-Ras fusion constructs which 
were treated with gold-labeled anti-GFP antibodies. The 
resulting immunogold point patterns were visualized with 
EM (immunoEM) to quantitatively describe Ras cluster- 
ing (Fig. [I]) 0- It was found that the classical raft 
model, wherein a fixed number of lipid rafts accommo- 
date a fixed fraction of raft-inserted proteins, is incom- 
patible with the observed gold point patterns [ll|. In- 
deed, for increasing expression levels, the classical model 
predicts an increase in the number of proteins per raft, 
and therefore a greater degree of clustering. To describe 
the data, Plowman et al. developed an alternative raft 
model, in which the size of Ras clusters remains constant. 
Assuming a constant, expression level- independent ratio 
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TABLE I: Representative Ras densities with corresponding 
numbers of Ras molecules on discretized lattice membrane as 
well as gold densities. Shown are the four Ras densities used 
in Figs. and [6] For lattice parameters, see Methods. 



of Ras in clusters to Ras monomers (cluster/monomer 
ratio) results in the formation of more rafts as expres- 
sion increases, and supports the notion that lipidated 
molecules such as Ras can drive the formation of rafts 
in order to create signaling platforms 
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native model predicts that 40% of active Ras molecules 
form clusters of radius 6-12nm, each containing about 
seven Ras molecules, and 60% are randomly distributed 
monomers While simulations of this model fit im- 
munoEM data, they do not provide a biophysical expla- 
nation for Ras clustering. Furthermore, these simulations 
violate laws of equilibrium physics. Specifically, the law 
of mass action predicts an increase in the fraction of clus- 
tered molecules as the exp ression level is increased (until 
membrane saturates) [ll[. This violation is troublesome 
as the experiments are done on in vitro membrane sheets, 
where no active, energy-driven processes can limit cluster 
size. Membrane sheets were fixed (and proteins immobi- 
lized) after membrane removal from cells Q, leading to 
equilibration of membrane and proteins prior to imaging. 

Recent experiments even go further and probe the de- 
sign principles of signaling by Ras clusters. Such studies 
suggest that, in the Ras-Raf-MEK-ERK pathway, Ras 
clusters act as an analog-digital- analog converter, where 
analog continuous EGF input is converted into digital, 
fully active clusters. The number of fully active Ras clus- 
ters, not the activity of individual Ras molecules, trans- 
lates into analog ERKpp output [2^, [2^. Specifically, 
these experiments show that Ras mutants with wide- 
ranging activities lead to the same total cellular ERKpp 
output [l^]. This suggests that Ras clusters act as digi- 
tal nanoswitches, which become fully activated even for 
small inputs. Furthermore, the concentrations of active 
Ras and ERKpp are directly proportional to EGF input 
[23 |. Hence, analog inputs produce analog outputs, me- 
diated by digital Ras clusters. 

Here we consider a physically-motivated model to 
study Ras clustering. The model mainly depends on a 
close-contact, attractive interaction between active Ras 
molecules (short range 2nm) and a repulsive interac- 
tion between Ras molecules irrespective of activity (long 
range ~ 5nm) . The short-range attraction promotes clus- 
tering of active Ras, while the long-range repulsion limits 
cluster size. Contrarily to previous models, we make no 
assumption about a constant, expression-level indepen- 
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FIG. 2: Relation between Lmax and Ras density A for im- 
munoEM data of gold labeled GFP-tH (black symbols) and 
RFP-tH (gray symbols), simulation averages and 99% confi- 
dence intervals (red), as well as a linear least-squares fit to 
simulation averages (red line). Lmax data points were ex- 
tracted from Ref. [ll[ with IMAGE J. (Left inset) gmax (black) 
and rmax (green) as a function of A. (Right inset) L max as a 
function of A without long-range repulsion (Vb = 0). Error 
bars represent standard deviations. For simulation details, 
including calculation of confidence intervals, see Methods. 



dent cluster/monomer ratio or the existence of lipid rafts, 
thus circumventing controversy surrounding their actual- 
ity. We equilibrate a discretized lattice membrane, occu- 
pied with active and inactive Ras molecules, using Monte 
Carlo simulations. After gold-labeling of Ras molecules 
from simulation outputs, wc perform a statistical cluster- 
ing analysis. The obtained statistical properties of Ras 
molecules quantitatively agree with the statistical prop- 
erties of immuno-gold point patterns for wide-ranging 
Ras expression levels (Fig. [2|) [ll]]. Our model makes 
predictions about the signaling properties of Ras clus- 
ters, supporting the notion that Ras clusters indeed act 
as an analog-digital-analog converter [23| . 



II. RESULTS 

Prior et al. studied Ras clustering in plasma 
membrane sheets using immunoEM of gold-labeled Ras 
molecules (Fig. [TIV). Gold point patterns were analyzed 
based on Ripley's K function. Specifically, the non-linear 
transformation L{r) — r was applied where r is the dis- 
tance between gold particles. This function is zero for 
complete randomness, positive for clustering, and nega- 
tive for depletion (Fig. [Tj3). Plowman et al. [ll| used 
the function's maximal value, termed Lmax for short, as 
summary statistics for clustering and found that Lmax is 
independent of Ras expression level (Fig. [51 symbols). 
This was rationalized by an ad hoc clustering model, as- 
suming a constant cluster/monomer ratio. Analysis of 
immuno-gold patterns is consistent with small clusters. 
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FIG. 3; Model ingredients. (A) Short-range attraction (red) 
and long-range repulsion (blue) as a function of distance be- 
tween two Ras molecules for the parameters given in Methods. 
Also shown is the cut-off beyond which the repulsive energy 
is set to zero (blue dashed line). (Inset) Representative part 
of lattice membrane showing three active Ras molecules (red) 
and one inactive Ras molecule (blue) . Neighboring active Ras 
molecules interact via the attractive short-range interaction 
(green bar). The cut-off used for the long-range repulsion 
is representatively shown for the central Ras (blue dashed 
circle). (B) Schematic of a gold-labeled antibody associated 
with a GFP-Ras molecule in the inner leafiet of the plasma 
membrane. 



containing approximately 6 to 8 molecules. Here we use a 
biophysical model of Ras clustering in the plasma mem- 
brane. In our model, a Ras molecule can be in either 
an active (on) or an inactive (off) state, corresponding 
to the respective GTP-bound and GDP-bound molecules 
for wild- type Ras. Both active and inactive Ras are as- 
sociated with membrane in line with experimental ob- 
servation [2^. The equilibrium probability of a single 
Ras molecule to be active depends on the effective free- 
energy difference between the on and off states, which in 
turn depends on input signals. We assume that active 
Ras molecules experience a short-range attraction, driv- 
ing cluster formation of active Ras, whereas a long-range 
repulsion limits cluster size (Fig. [5J\, main panel). Such 
a long-range interaction may result from lipid-anchor in- 
duced membrane deformations. To obtain equilibrium 
properties, we approximate the membrane by a square 
lattice, populated by Ras molecules of a specified den- 
sity (Fig. inset), and perform Monte Carlo simu- 
lations. For comparison with immunoEM experiments, 
we added lOnm-long gold-labeled antibodies (maximally 
one per Ras) to the Ras molecules in the experimentally 
observed capture ratio (Fig. [2J3). We mainly use the four 
Ras densities given in Table |T1 To specifically compare 
with experiments on varying Ras-expression level (sym- 
bols in Fig. [21 main panel), we calculate the Lmax value 
for additional Ras densities. Note that these experiments 
are based on the lipid anchor of H-Ras (tH), which has 
similar clustering properties as active H-Ras [11| . For de- 
tails on the experiments and our approach, see Methods. 

Figure [¥] shows typical equilibrated membrane lattices 
for the four Ras densities (left panels) with the corre- 
sponding plots of L(r) — r (right panels). For the low- 
est density, individual L(r) — r plots are highly variable. 
To produce meaningful statements about clustering we 
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FIG. 4: Monte Carlo simulations and point-pattern analysis. 
Snapshots of equilibrated Ras molecules on lattice membrane 
(left column; active Ras in red and inactive Ras in blue) and 
corresponding L{r)—r plots (right column) after gold labeling 
for the four densities from Table |T] (density of Ras molecules 
increases from top to bottom). Shown in the L{r) — r plots 
are individual simulations (cyan curves) , their averages (thick 
black curves), as well as 68.3%, 95.4%, and 99.0% confidence 
intervals (red, green, and blue dashed lines, respectively). For 
simulation details, including calculation of confidence inter- 
vals, see Methods. 



also show the averaged plot, as well as provide confi- 
dence intervals. In line with experiment on varying Ras- 
expression level, we observe that for our model, L„iax is 
approximately independent of Ras density (Fig. [21 main 
panel) . The same is true if the analysis is done directly on 
Ras molecules instead of the gold particles, demonstrat- 
ing the robustness of the result with respect to the details 
of Ras labeling by gold. Distance rmax, defined as the dis- 
tance corresponding to L„iax, is about 8nm (Fig. [21 left 
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FIG. 5: Monte Carlo simulations and point-pattern analysis 
for conventional clustering model without long-range repul- 
sion (Vb = 0). For a description of symbols and lines, see Fig. 
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inset), or equivalently, 4 Ras molecules. Hence, clusters 
contain few, about 4 to 10, Ras molecules. An alternative 
clustering analysis based on the pair-correlation function 
g{r) gives similar results, i.e. gmax values are indepen- 
dent of A (Fig. [21 left inset). Hence, cluster sizes and 
their dependence on expression level are in good agree- 
ment with previous estimates [lTj |. 

We also explored a more conventional clustering model 
without the long-range repulsion, but maintaining the 
short-range attraction. As shown in Fig. [5] (left pan- 
els), Ras molecules form increasingly larger clusters at 
increasing Ras densities. Examination of the L{r) — r 
plots (right panels) shows that imax decreases for in- 
creasing Ras densities (Fig. [21 right inset), which is in 
stark contrast to experiments. Hence, limiting the cluster 
size by the long-range repulsion is a necessary ingredient 
to correctly describe immunoEM data and, hence, Ras 
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clustering. 

Next we examined the fraction of Ras molecules in clus- 
ters. Previous models assumed that the fraction is con- 
stant, i.e. independent of Ras density. In contrast, Fig. 
ini shows that in our model the distribution of the fraction 
clustered increases significantly with density, indicating 
that a constant cluster/monomer ratio is not required to 
describe the immunoEM data in Fig. [2] Also shown in 
the Fig. [5] is the fraction clustered for the conventional 
clustering model without the long-range repulsion. The 
distribution also shifts to higher values with density, al- 
though to a lesser extent as fractions arc much higher 
to start out with due to the missing long-range range 
repulsion. To clearly rule out the conventional cluster- 
ing model as a suitable model, we tested whether as- 
suming a constant fraction clustered (or equivalently, a 
constant cluster/monomer ratio) can explain the immu- 
noEM data. For this purpose we collected simulations 
from different densities but same fraction clustered and 
compared their imax values. However, even with this 
strong selectivity of simulations, imax values continued 
to decrease with increasing density (Fig. [6l inset). 

Figure[7]shows the signaling characteristics of Ras clus- 
ters for four different inputs. For input we use the free- 
energy difference between on (active) and off (inactive) 
Ras states (c/. Eq. [T]). To test if our model produces 
digital-like nanoswitches, which are fully active even for 
small inputs, we identified clusters of two or more con- 
nected Ras molecules and calculated the cluster activ- 
ity, i.e. the fraction of active Ras molecules in clusters. 
The bar chart in Fig. [TlA. shows that our model indeed 
produces nanoswitches, which are fully active even for 
small stimuli, as indicated by experiments (27| . In con- 
trast, the activity of a single Ras molecule does not be- 
have like a switch (Fig. [8lA., dashed line). Furthermore, 
Fig. [7J3 provides the total activity of all Ras molecules 
in the membrane irrespective of whether Ras molecules 
belong to clusters or not. We find the total activity is 
approximately proportional to the input (black line) in 
the range considered here in line with experiment [2J]. 
In our model, this is due to the fact that the number of 
Ras clusters is proportional to the input (Fig. [7j3, blue 
line). 

There has recently been immense interest in under- 
standing the effect of noise in signal transduction [29l.[30t. 
Biochemical reactions are inherently noisy as they are 
based on random collisions of molecules. This intrin- 
sic noise is further enhanced by the small number of 
molecules involved. Furthermore, rate constants may 
fluctuate, as they depend on external conditions such 
as other molecules not explicitly considered as part of 
the biochemical reactions. This extrinsic noise also in- 
cludes fluctuations in the input itself. To address how 
Ras signaling is affected by noise, we compare signaling 
by Ras clusters (Fig. [7|) with signaling by non-interacting 
Ras molecules without clustering ability (Fig. [8|). In- 
trinsic noise is inherently part of our simulations as Ras 
is allowed to randomly switch between the active and 
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FIG. 6: Distributions of Ras fractions in clusters. Different 
colors correspond to the four Ras densities from Table IT] i.e. 
A = 625 (red), A = 1250 (blue), A = 2500 (green), A = 5000 
(black) in units of /im~^. Shown are results with (solid lines) 
and without (dashed lines) long-range repulsion. A Ras clus- 
ter is defined as two or more connected Ras molecules. (Inset) 
imax for pairwise constant fractions (overlapping fractions), 
i.e. fraction range 0.72-0.75 for A = 625 and 1250 (circles and 
dashed line), fraction range 0.81-0.84 for A = 1250 and 2500 
(triangles up and dotted line), and fraction range 0.88-0.91 for 
A = 2500 and 5000 (triangles down and dashed-dotted line). 



the inactive states. The intrinsic noise for the activity 
of Ras in clusters (Fig. [Tl'V, black error bars) is signifi- 
cantly less than for the activity of a single Ras molecule 
(Fig. black error bars) since clusters are fully active 
and hence suppress random switching. This difference in 
intrinsic noise is reduced when considering the intrinsic 
noise of the total activity from all Ras molecules in the 
membrane, which is only slightly smaller for Ras clus- 
ters (Fig. [7J3, black error bars) than for non-interacting 
(Vo = J = 0) Ras molecules (Fig. [5J3, black error bars). 
This is due to the fact that the number of Ras clus- 
ters, which is necessarily smaller than the number of Ras 
molecules, can fluctuate significantly (Fig. [7(3, blue error 
bars). Most importantly, Ras clusters are more robust to 
input noise, at least for sufficiently large inputs, than 
non-interacting Ras molecules (by comparison of green 
error bars in Fig. [TjA. and Fig. [H1A-). Here, input noise 
represents fluctuations in input much faster than assem- 
bly/disassembly of clusters but slower than Ras signal- 
ing. Therefore, it is assumed that extrinsic noise only 
affects the activity of Ras molecules, not clustering itself 
(see captions of Fig. [TjA. and Fig. [SIA. for details). This 
shows that Ras clusters have superior signaling proper- 
ties compared to non-interacting Ras molecules without 
clustering ability. 
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FIG. 7: Signaling properties of Ras clusters. (A) Cluster 
activity as a function of input (parameter Ae). Cluster activ- 
ity is defined as fraction of active Ras in clusters from sim- 
ulations (bar chart), where a cluster contains two or more 
contacting Ras molecules. Also shown is approximate clus- 
ter activity Pon = 1/[1 -|- exp(7VAe)], which assumes that 
all A'^ Ras molecules in a cluster (here chose N = 10) are 
tightly coupled and hence are either all on (active) or all off 
(inactive) together (black line). Black error bars show stan- 
dard deviation and represent intrinsic noise. Green error bars 
represent extrinsic noise, calculated with noise propagation 
formula SPon = (dPon/dAe)5[Ae] = APon(l - Pon)5[Ae] for 
S[Ae] = O.OSfcsT. (B) Total activity of aU Ras molecules 
in the membrane, normalized by the total number of Ras 
molecules (grey bar chart, left axis), and number of Ras clus- 
ters (white bar chart, right axis). Also shown are linear fits. 
Error bars represent standard deviations. To enlarge black 
error bars for better visualization in B, we used the square- 
root of the total variance from pooled simulations of inputs 
Ae, Ae -I- O.lksT, and Ae - O.lfcsT. 



III. DISCUSSION 

Different Ras isoforms are known to form nonoverlap- 
ping signaling clusters [1, important for localized 
signaling of the Ras-Raf-MEK-ERK pathway [HI, HI, 
involved in cell proliferation, differentiation, and apop- 
tosis [20| . In addition to the fundamental importance of 
Ras in this pathway, Ras mutations are found in 30% of 
human cancers [l^, [3, [3] ■ Ras clusters are also consid- 
ered evidence of lipid rafts [lH . Lipid rafts have attracted 
considerable interest due to their alleged role in protein 
sorting and specificity of signaling [l|, 0, 3 ■ this work, 
we provided a biophysical model of Ras clustering, and 
compared results with gold-point patterns obtained from 
immunoEM of plasma membrane extracts (Fig. [T|). In 
particular, we obtained that clustering of Ras molecules, 
i.e. the cluster/monomer ratio, is independent of expres- 
sion level (Fig. [2]), in line with experiments on the lipid 
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FIG. 8: Signaling properties of non-interacting Ras molecules. 
(A) Activity of single Ras molecule (dashed line; calculated 
with Eq. [T]for Pon) as a function of input (parameter Ae). 
Black error bars represent intrinsic noise, calculated from the 
square-root of the binomial variance Pon(l — Pon). Green er- 
ror bars are approximately 0.01 in magnitude and represent 
extrinsic noise, calculated with the noise propagation formula 
SPon = Pon(l - Pon)(5[Ae] aud 5[Ae] = O.OSfcsT. (B) Total 
activity of all Ras molecules in the membrane, normalized by 
the total number of Ras molecules (bar chart) and linear fit 
(dashed line) . Error bars represent standard deviation, calcu- 
lated from the square-root of the total variance from pooled 
simulations of inputs Ae, Ae + O.lfcsT, and Ae — O.lfcsT. 



anchor of H-Ras (tH) [Tl| . In our model, as well as in ex- 
periments, clustering is quantified by the maximum value 
(termed Lmax) of function L{r) — r where r is the 
distance between gold particles. Our model has two main 
ingredients exemplified in Fig. [31 (1) a short-range at- 
traction between active Ras molecules (e.g. Ras-GTP) 
promoting clustering, and (2) a long-range repulsion be- 
tween Ras molecules, which limits cluster size. 

Another important feature, which makes our model 
fundamentally different from previous Ras clustering 
models [HI [31, is that the fraction of clustered Ras 
molecules is not a model parameter [llj but a prediction 
from our simulations. Indeed, if we calculate the fraction 
of clustered molecules for the four densities from Table 
m we obtain the distributions shown in Fig. [6l The frac- 
tion of clustered Ras increases with density, indicating 
that the assumption of a constant cluster/monomer ratio 
[ill . [2^ is misleading for descri bing the immunoEM data 
for different expression levels [llj . Since this assump- 
tion violates equilibrium thermodynamics, our model is 
more suitable for describing in vitro immunoEM data in 
absence of energy sources from the cell. Note that in 
living cells clustering may in part be regulated by ac- 
tive, energy-dependent mechanisms. For instance, clus- 
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tering and signaling of constitutively active K-RasG12V 
depends on the presence of actin fences Such 
membrane-associated actin filaments, part of the actin 
cortex, are highly dynamic and, hence, K-Ras cluster- 
ing may be regulated. An expression Icvcl-indcpcndcnt 
cluster/monomer ratio has also been found for glycosyl- 
phosphatidylinositol-anchored proteins (GPI-AP) in vivo 
[33|. Finally, clustering of proteins in the immunological 
synapse is an active, actin-myosin dependent process [3J] , 
presumably to overcome the entropic barrier of localizing 
proteins [35[- 

What is the role of Ras clusters beyond simple pro- 
tein sorting? Harding et al. argued that Ras clusters 
allow for highly precise coding of time-dependent in- 
puts, termed high fidelity signaling |25i, i26||. First, Ras 
is highly abundant in the membrane (tens of thousands 
molecules), hence the number of active clusters can be 
wide-ranging depending on input, e.g., EGF. Second, 
clusters have a short life time (about 0.4 s). This high 
turnover of clusters allows for high temporal precision. 
Third, Ras clusters act as digital nanoswitches, which 
may lead to noise reduction in the signal transmission 
step across the membrane due to coarse graining and av- 
eraging of rapidly fluctuating signals. In fact, Ras clus- 
ters have similarity to an analog-digital-analog converter 
known from engineering, transmitting analog EGF input 
into analog ERKpp output [2^, [2^. This design prin- 
ciple was indeed recently confirmed by experiments on 
the Ras-Raf-MEK-ERK pathway in baby hamster kidney 
cells. In particular, experiments showed that Ras clus- 
ters (actually Raf*-tII with varying kinase activity) are 
fully active even for small inputs (kinase activities) [131, 
and that the total amount of active Ras (concentration 
of K-Ras-GTP) and ERKpp are proportional to EGF in- 
put [2^. Such analog ERK activation was recently also 
observed in proliferating mammalian fibroblasts (36[ . 

The above listed properties of Ras clusters are sup- 
ported by our model. According to Fig. Ras clus- 
ters arc fully active, even for small inputs. Nevertheless, 
the number of Ras clusters and hence the total activ- 
ity of all Ras molecules in the membrane are approx- 
imately proportional to the input (Fig. [7J3), allowing 
faithful transmission of continuous, time-dependent in- 
put signals. Interestingly, the activity of a single Ras 
molecule and of non-interacting Ras molecules are also 
approximately proportional to the input (Fig. [H]A and 
Fig. [5J3, respectively). However, signaling by Ras clus- 
ters is less noisy and, hence, Ras clusters can transmit 
signals more robustly than non-interacting Ras molecules 
without clustering ability. The activity of Ras in clusters 
exhibits smaller intrinsic noise from random switching 
between active and inactive states (Fig. [7l?V and B, black 
error bars). The activity of Ras in clusters is also less 
sensitive to extrinsic noise from fast fluctuations in in- 
put than non-interacting Ras molecules (by comparison 
of green error bars in Fig. [Tl'V and Fig. respec- 
tively). The reason for the noise reduction by clusters 
is that Ras clusters are fully active, suppressing random 



switching between active and inactive states, as well as 
activity changes due to fluctuations in input. 

Our model relies on the short-range attraction and 
long-range repulsion between Ras molecules. The physi- 
cal origin of these interactions are yet to be determined. 
However, the attraction may originate from direct Ras- 
Ras interaction via hy drop hobic, van der Waals, or elec- 
trostatic interactions [37|, but may also be mediated 
indirectly by scaffold proteins and lipids [13, [H, [37| . 
The latter mechanism is supported by the finding that 
the positively-charged polybasic C-terminus of K-Ras 
binds negatively charged phospholipids and sequesters 
acidic phospholipids, which may attract even more K-Ras 
molecules. Furthermore, that mutant GFP-K-RasG12V 
S181D has a reduced ability to bind the membrane, as 
well as to cluster |37| . Such a lipid- mediated mechanism 



would support the concept of dynamic lipid rafts, which 
only form in presence of an activated membrane protein 
such as Ras-GTP. The physical origin of the repulsion is 
harder to pinpoint, but may result from induced mem- 
brane curvature as a result of insertion of the farnesyl- 
polybasic anchor of K-Ras or the farnesyl-palmitate an- 
chor of H-Ras into the inner leaflet of the membrane ■ 
Lipid-anchor induced membrane deformations are sup- 
ported by molecular dynamics simulations [s^ , and may 
lead to long-range repulsion [s^ . Furthermore, recent 
experiments explicitly show that smallGTPases (Arf ) in- 
duce membrane curvature (40| . 

There are certain shortcomings of the immunoEM 
data, rendering the cluster analysis of immuno-gold data 
difficult. In our simulations, the addition of gold parti- 
cles to completely random distributions of Ras molecules 
produced the following. If more than one gold-labeled an- 
tibody is allowed to bind a Ras molecule, provided there 
are no steric clashes between gold particles, the L{r) ~ r 
plot still predicts Ras clustering. Varying the antibody 
length systematically resulted in distance Tmax being ap- 
proximately equal to the length of the antibody. This is 
due to a small fraction of Ras molecules being associated 
with multiple gold particles: these gold particles will be 
within two antibody lengths of each other, and approxi- 
mately one antibody length from each other on average, 
resulting in a peak in the L(r) — r plot. This suggests 
that, unless it is verified that each Ras can only bind 
a single antibody (i.e. the anti-GFP antibody can only 
bind a single epitope on GFP fused to a Ras molecule), 
immunoEM data can overestimate clustering. In con- 
trast, two clustered Ras molecules in contact with each 
other could each interact with separate antigen-binding 
regions of the same antibody, since an antibody has two 
antigen-binding regions. In this case the cluster of two 
Ras molecules is unobservable by immunoEM, leading to 
an understimation of clustering in the gold point pat- 
terns. These issues would have to be addressed if im- 
munoEM studies are to form the basis of an accurate 
quantification of Ras clustering. 

In conclusion, a comprehensive description of Ras clus- 
tering is an essential step in the understanding of Ras 
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signaling properties, and of small, inner-membrane GT- 
Pases in general. For instance, we showed that clustering 
leads to robustness to noise, especially input noise (Figs. 
[7] and While the model we have analyzed fits the 
data (Figs. [2] and [7]), several questions remain unan- 
swered. For one. lack of high resolution structural infor- 
mation about clustered Ras molecules prevents us from 
"seeing" clearly into the physicochemical basis of clus- 
tering (only partial crystal structures of Ras molecules 
exist l41[). There is also the possibility of regulation of 
clusters from within utilizing specific lipids and scaffold 
proteins, which has scarcely been addressed in the litera- 
ture thus far [l^ , but would provide critical details to the 
construction of an accurate model for clustering. While 
we have shown that immunoEM data can aid in the vi- 
sualization of Ras clusters, data regarding the dynamics 
of_clustering are found in the form of SPT FRET 
fluorescence recovery after photobleaching (FRAP) 
and single-molecule fluorescence microscopy stud- 
ies These data remain to be integrated into 
more detailed spatio-temporal Monte Carlo simulations, 
so that the exchange of proteins between freely diffus- 
ing monomers in the membrane and immobile clusters 
can be investigated [ij. Interestingly, our biophysical 
model of Ras clustering shares the long-range repulsion 
due to elasticity with recent models of lipid microphase 
separation [3^ and chemoreccptor clustering (45l . [4^ in 
bacteria. Hence, similar biophysical principles may gov- 
ern the clustering of very different types of proteins in 
prokaryotic and eukaryotic membranes. The latter may 
include EGF and Fc7 receptors, which are believed to 
associate with rafts or to form small clusters StI HM . 



IV. METHODS 

Experimental immunoEM data 

Relevant experiments are described in 0, [O, [13, US- 
Briefly, Ras clustering was examined on intact 2-D 
sheets of apical plasma membrane, ripped off from 
adherent baby hamster kidney cells directly onto EM 
grids. Ras-fluoresccnt protein fusion constructs were 
used, including the minimal plasma membrane tar- 
geting motif (lipid anchor) of H-Ras fused to GFP 
(GFP-tH) or RFP (RFP-tH), as weU as constitutively 
active H-RasG12V and K-RasG12V. These are tagged 
using affinity-purified polyclonal anti-GFP antibodies, 
conjugated with 4nm gold particles, and visualized using 
electron microscopy (immunoEM). The resulting point 
patterns of gold particles were analysed for clustering 
(see below for details) . Ras isoform clustering was found 
to depend differentially on membrane-associated actin 
lipid-raft constituent cholesterol Q, and scaffold 
proteins galectin-1, galectin-3, and Sur-8 [i^. For a 
recent review see UM ■ 



Biophysical model 

A Ras molecule in the membrane can be in either one 
of two states, active (on) with energy eon or inactive 
(off) with energy [50|. For wild-type Ras, the active 
(inactive) state corresponds to Ras-GTP (Ras-GDP). 
More generally, the two states correspond to two differ- 
ent protein conformations, making the two-state model 
applicable to activity mutants and lipid anchors as well. 
For any such two level system, the probability for a 
single Ras molecule to be active is 



p — 

^ on 



1 



1 



(1) 



where Ae = Con — eoff is the free-energy difference between 
the active and inactive states. While Eq. [l]is not explic- 
itly used for our simulations as it describes the activity of 
Ras in absense of interactions, it builds intuition about 
parameter Ae. This parameter is effectively determined 
by the input signal of the pathway, e.g. EGF, except 
for the Ras activity mutants and lipid anchors, where it 
desribes an energetic bias in conformational state. 

To describe clustering of active Ras molecules, as di- 
rectly observed for K-Ras and H-Ras using in vivo FRET 
[2]| . we introduce short-range attraction J between ac- 
tive Ras molecules, driving cluster formation. In order 
to limit cluster size, we introduce long-range repulsion 
V{r), where r is the distance between two Ras molecules. 
For the repulsive interaction energy, we use a Gaussian 
function as previously applied for describing microphase 
separation of lipid mixtures [s^ 



yif) = Vq ■ exp 



2a2 



(2) 



where Vq is the maximal repulsion for two Ras molecules 
in close proximity and a is the width, i.e. the range 
of the repulsion beyond which the potential drops 
quickly (Fig. [SK). The frustration between short- 
range attraction and long-range repulsion leads to 
small clusters. More precisely, the optimal cluster 
size corresponds to the minimum of the cluster energy 
divided by the number of Ras molecules in the cluster, 
i.e. the energy density [5l] |. The parameters used 
in this study are Ae = -0.8, J = -5.0, Vq = 2.0, 
and a = 2nm. Long-range repulsions were neglected 
beyond a 6nm cut-off to reduce the calculation time. All 
energies are in thermal energy units ksT with ks being 
the Boltzmann constant and T the absolute temperature. 

Monte Carlo simulations 

Since the immunoEM data are obtained from in vitro 
membrane sheets, clustering is an equilibrium process. 
Models of such phenomena are therefore particularly 
amenable to Monte Carlo simulations, which include en- 
ergetics as well as entropy [5lj. To set up simulations, we 
discretize a defined area of the plasma membrane inner 
leafiet to obtain a two-dimensional M x M square-lattice 
where M is the lattice size. On the lattice, each position 
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is uniquely described by an index i (if the 2-D lattice is 
thought of as a linear array of length M^). We assign a 
Boolean value Si to every Ras, where = if Ras i is 
active, and = 1 if Ras i is inactive. Using this nota- 
tion, we can construct an energy function describing the 
total energy E for a set of N molecules 

iV 

i? = ^Ae.(l-s,) + ^ ./•(!- s,)(l-s,) + ^y(r,,) 



where {i,j) denotes nearest-neighbor pairs. 



(3) 



After randomly generating the positions of the starting 
Ras molecules on the lattice, individual Ras molecules 
are chosen at random and attempted to move to a 
new location on the lattice. Included in each step is 
a probability of switching between active and inactive 
Ras. Moves are accepted or rejected based on the 
Metropolis-Hastings algorithm. We use a lattice of size 
M = 300 and a lattice constant a = 2nm (the size of a 
Ras molecule [31), resulting in a 0.36 ^m^ membrane. 
In order to reduce boundary effects of the lattice, we 
adopt periodic boundary conditions. 



Gold particles 

In order to compare the simulation outputs with immu- 
noEM data, gold-labeled antibodies are added to the 
equilibrated Ras molecules. The length of the antibody 
used in the experiments is lOnm [llj. When an antibody 
binds a Ras molecule, the gold particle associated with 
the antibody can at any one time occupy any position 
on the surface of a hemisphere around the Ras molecule 
(Fig. [3j3). For simplicity, the radius of the hemisphere 
is chosen equal to the length of the antibody, and 
the centre of the gold particle's position is projected 
onto the plane of the membrane, defining the particle's 
position on the lattice. This position is then matched 
against previous gold positions for steric clashes, and 
if it is found to be closer than 4nm from another gold 
particle, the position is rejected and another Ras is 
picked at random. This process is iterated until 42% 
of Ras molecules are occupied, corresponding to the 
experimentally-observed capture ratio 



Cluster analysis 

We use three functions to evaluate the degree of Ras clus- 
tering: K{r), L{r)~r and g{r). Ripley's -fsT-function iC(r) 
was first proposed for analyzing spatial point patterns 
[5^ . K{r) calculates the expected number of particles 
within a distance r of any particle, normalized by the 



average density A 
N{r) 



K(r) 



A 
1 1 



N N 



(4) 
(5) 



1=1 i^j 



A 



N N 



N{N 



-^T.J2^^^i\\^,-^^\\<r) (6) 



1=1 i^j 



where A is the area of the lattice studied, N the number 
of Ras molecules or gold particles, A the surface density 
and lij (x) an indicator function which takes a value of 1 if 
\\xj —Xi\\ < r and otherwise. Under the null hypothesis 
of complete spatial randomness, N{r) ~ Avrr^, so K(r) = 
7rr^. An often used non-linear transformation of K(r) 
which we shall employ is 0] 



L(r) — r = 



K{r) 



(7) 



which has a value of for complete spatial randomness, 
is positive for clustering, and negative for depletion of 
particles. For large r, L{r) — r is zero on average, since 
particles are uncorrelated. Since L{r) — r is a non-linear 
transformation of K{r), when averaging over multiple 
simulations to resemble a large piece of membrane (see 
below), the K{r) values are averaged first and only then 
transformed into L{r) — r. For further analyses of simu- 
lations, we use summary statistics Lmax = max[L(r) — r] 
and corresponding distance r^ax = argmaxr[L(r) — r]. 

The pair-correlation function g{r) can be defined in 
two different ways. The first by normalizing and differ- 
entiating K{r) [5I, HI] 



1 dK{r) 
Itxt dr 



(8) 



and the second by counting in a similar manner to Kir) 
but in concentric rings: 

A ^ N N 

(9) 

for r > a, where a is the lattice constant. Testing 
these two versions of the pair-correlation function yielded 
slightly different absolute values of g{r), but the relative 
behaviors of the two were identical. For a random distri- 
bution of particles and for large r in general, g(r) takes a 
value of 1 on average. We again use summary statistics 
ffmax = max5(r) and r^ax = argmax^ .g(r). 

To estimate confidence intervals for the L(r) — r 
cluster analysis, 99 simulations were run for each density 
with all interactions set to zero, simulating a random 
distribution of Ras in order to obtain an estimate of the 
background clustering noise intrinsic to each density. 
Triplets of K{r) were averaged to simulate a 1/im^ 
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membrane as used in experiments, and L{r) — r values 
were calculated for each. The 68.3%, 95.4%, and 99.0% 
confidence intervals for individual measurements were 
obtained by multiplying the standard deviation of the 
33 triplets by 1, 2, and 2.576, respectively. Standard 



deviations for L,„ax, ffmax, 
based on triplets as well. 



and Tn 



were calculated 



the proteins mentioned in the text are: H-Ras P01112, 
K-Ras P01116, and N-Ras POllll. 
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